Excitonic Dynamical Franz-Keldysh Effect 
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The Dynamical Franz-Keldysh Effect is exposed by exploring near-bandgap absorption in the pres- 
ence of intense THz electric fields. It bridges the gap between the DC Franz- Keldysh effect and 
multi-photon absorption and competes with the THz AC Stark Effect in shifting the energy of the 
excitonic resonance. A theoretical model which includes the strong THz field non-perturbatively 
via a non-equilibrium Green Functions technique is able to describe the Dynamical Franz-Keldysh 
Effect in the presence of excitonic absorption. 
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We report observation of the Dynamical Franz Keldysh 
Effect (DFKE) by examination of the near-bandgap opti- 
cal properties of a semiconductor multiple quantum well 
(MQW) in a new experimental regime. By applying a 
strong electric field at frequencies near 1 THz, we per- 
form electro-absorption spectroscopy which is described 
by neither the DC Franz-Keldysh Effect (FKE) @J|, nor 
by optical-frequency effects like multi-photon absorption 
(MPA) ||. The DFKE region, bridging these two ex- 
tremes, is a topic of longstanding theoretical interest 
but has not previously been studied experimen- 
tally. 

The crossover between "high" and "low" frequency 
may be quantified in terms of two parameters, .Eke an( i 
7, defined as: 
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Eke is the mean kinetic energy of a particle of mass 
to* and charge e in an electric field -Ethz cos(wxHzi), i-e- 
the ponderomotive energy. The DC FKE corresponds 
to 7 ^> 1; multi-photon effects occur for 7 <C 1; and the 
case of 7 ~ 1 is a new, distinct DFKE regime Q (Fig. 1). 
This regime is difficult to access experimentally, as Eke 
must be comparable to or greater than any broadening 
energies in the system. However, it is possible for applied 
fields with w T Hz ~ ITHz and E THz - 1 - lOkV/cm §. 

In this paper, we explore the optical absorption of a 
semiconductor multiple quantum well (MQW) driven by 
an intense THz electric field polarized in the plane of the 
MQW layers. For the case of non-interacting carriers, 



the DFKE predicts a blue-shift of the main absorption 
edge in the spectrum by Eke, as well as increased sub- 
gap absorption §] (Fig. 2a). This DFKE blue-shift has 
not previously been observed. 

In order to properly address the DFKE, we must also 
consider the effect of excitons, which dominate the ob- 
served near-bandgap absorption in undoped semiconduc- 
tors. Excitons are bound states of an electron and hole 
with hydrogen-like quantized energy levels. The energy 
level spacings can be in the THz regime. The interaction 
of excitonic states with a THz field leads to shifts of the 
levels via the AC Stark Effect § (Fig. 2b). The THz field 
interacts most strongly with the exciton Is — ■> 2p tran- 
sition j8|. We present below a novel theoretical model 
which shows that the experimentally observed changes 
in the absorption spectrum, as a function of the THz 
field intensity and frequency, display the DFKE and its 
interplay with the AC Stark Effect. 

Our MBE-grown sample consists of a 20-period MQW 
with 80A In . 2 Gao.8As wells separated by 150A GaAs 
barriers on a GaAs substrate. The energy of the 2D 
MQW states is below the substrate band-gap, allowing us 
to perform direct transmission spectroscopy without sub- 
strate removal. Strain lifts the valence-band degeneracy 
at k = 0, splitting the light-hole and heavy-hole bands 
by 50meV. The electron-heavy-hole exciton is the lowest- 
energy feature. It is inhomogeneously broadened by in- 
terface roughness and alloy disorder || with a FWHM 
of 6meV. The exciton Is — > 2p transition energy (/KJ12) 
is ~8meV, estimated from magneto-transmission experi- 
ments ||. 

A CW Ti-Sapphire laser provides tunable near- 
bandgap radiation, in the near infrared (NIR), focussed 
to a ~ 200^m spot. Its transmission is detected by a con- 
ventional Si PIN photo-diode, and peak NIR intensity is 
~ 100mW/cm 2 . THz radiation is provided by the UCSB 
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Free-Electron Laser, a source of intense, coherent radia- 
tion tunable from ftwxHz ~0.5-20meV It is roughly 
collinear with the NIR beam and focussed to a 0.5-2. 5mm 
spot coincident with the NIR focus (Fig. 2c). We mea- 
sure transmission spectra as a function of THz intensity 
and frequency. All measurements are performed in He 
vapor at ~7K. 

We observe two distinct regimes of experimental be- 
havior for different ranges of THz frequency, wthz ■ When 
wthz is less than W12, we observe a red-shift of the exci- 
ton for low THz intensities, i.e., an AC Stark shift. As 
THz intensity increases, this shift reaches a maximum 
and reverses, eventually becoming a net blue-shift at the 
highest THz intensity; the DFKE blue-shift dominates 
at high intensities as the AC Stark Effect saturates (Fig. 
3a). Conversely, if wthz is greater than tui2, we observe 
a blue-shift of the exciton which increases monotonically 
with increasing THz intensity; the AC Stark Effect and 
DFKE act in concert, each contributing to the blue-shift 
(Fig. 3b). In both cases, the exciton peak is broadened 
and suppressed, and the broadening increases with in- 
creased THz intensity (Fig. 3a, b). 

Since the apparent peak shifts are small, we quantify 
them by taking a weighted average of the data near the 
exciton line. We find the geometric mean of the upper 
half of the excitonic peak, and plot this "center of mass" 
against 7, which is directly proportional to THz power 
(Fig. 4). Note that our estimates of the maximum 7 at 
each frequency are only accurate to within a factor of 2 
or 3. We also show results of our theoretical calculations 
(to be discussed below). Again, we see an initial red- 
shift, followed by a blue-shift at low wthz (Fig. 4a) and 
a monotonic blue-shift at high wthz (Fig. 4b). 

Theoretically, we describe the optical properties of the 
system by calculating the time dependent induced polar- 
ization. The system is probed using a weak NIR probe 
field -Enir(£) with frequency u> ~ E sap /h. To linear or- 
der in the NIR field, the macroscopic polarization is the 
/c-trace of 
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All the information about the THz field is contained in 
the non-equilibrium generalization of the retarded two- 
time interband susceptibility x r (fc;i,i'). In our previ- 
ous work based on non-equilibrium Green functions, 
we studied the susceptibility starting from a two-band 
Hamiltonian, including the THz field non-perturbatively 
in the Coulomb-gauge via the vector potential A(t) = 

— Ethz sin(ciJTHzi)/wTHz- Here, -Ethz is oriented in the 
plane of the MQW layers. Non-equilibrium Green func- 
tion techniques Jljj yield an analytic expression for the 
single particle susceptibility: 
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where d is a dipole matrix-element and e[k] = ^ -r j^gap, 

fi being the electron-hole reduced mass and k the in-plane 
momentum. 

We include excitonic effects by using the Bethe- 
Salpeter equation to determine a susceptibility which 
includes the electron-hole interaction in the ladder ap- 
proximation: 

x r (k;t,t')=x r (k;t,t') 

d 2 P r°° - - 

dt" x r (k;t,t")V(\k-k'\) X r (k';t'\t'), (5) 
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where x r (k',t, t') is the single-particle, non-interacting 
susceptibility of Eq. (||). Expressing the susceptibility 

as X r (k;t,t') = E„/Jf xUk^)e^ t - t '^ n2 ^ t + t '\ 
the integral equation (0) becomes a matrix equation [0| , 



Xn(k,Uj) = Xn(k,U>) 

+ ^2Xn-ni{k,U + 2n'uJ THz ) 
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-V(\k-k'\) X r n ,(k , ,LJ + 2(n' -n)w T Hz). 



For physically achievable THz intensities, only small val- 
ues of n are needed. 

Numerically, we solve the matrix equation including 
both s-wave and p-wave scattering; both are important, 
since the THz field will couple the s and p states of the 
exciton and influence the observed resonance via the AC 
Stark Effect. The solutions to the resulting equation are 
found by discretizing the integrals to yield a set of linear 
equations which we solve by standard methods |14|. 

Finally, we relate the results to physically measurable 
quantities by expressing the macroscopic polarization as 



P(t) = E Nm J2xn{u)e 
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where Xn(^>) = J2kXn(k,^)- The linear absorption is 
thus proportional to Im%o(<^)- The terms with n ^ de- 
scribe the non-linear mixing of the NIR and the THz field, 
resulting in optical sideband generation. These sidebands 
show a rich structure, but an analysis is beyond the scope 
of this work, and will be reported separately. 

Using the above method, we have calculated absorp- 
tion spectra using the experimental parameters of our 
system. We find good agreement with experiment, us- 
ing no fitted parameters. For wthz < ^12, theory pre- 
dicts a red-shift of the exciton absorption peak at low 
THz intensities. With increasing THz intensity, the shift 
saturates and then reverses, eventually becoming a net 
blue-shift (Fig. 5a). Theory also predicts a blue-shift for 
wthz > ^12) which monotonically increases with THz in- 
tensity for experimentally- accessible intensities (Fig. 5b). 
In both cases, the exciton peak is suppressed and inhomo- 
geneously broadened with increasing intensity, i.e., with 
increasing 7 (Fig. 5ab). Using the same method as for 
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the measured spectra, we compute the "center of mass" 
of the calculated exciton lines. We plot the shift of the 
theoretical "center of mass" as a function of 7 in Fig. 
4, along with the experimental results. The theoretical 
and experimental peak shifts show identical qualitative 
behavior. 

As described above, this behavior can be understood 
as the DFKE acting in competition with the AC Stark 
Effect. The AC Stark Effect, at low -Ethz > results in a 
shift of the exciton level, 
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where V is the width of the Is — > 2p transition line [Q. 
We estimate fiF ~ 4meV in our sample. The magnitude 
of A is dominated by T, and we expect little resonant 
enhancement of |A| for wthz ~ ^>12- A does not in- 
crease indefinitely, and will saturate as -Ethz increases, 
going from quadratic to linear dependence on £thz 0- 
The DFKE, however, always provides a blue-shift pro- 
portional to -Ethz ! given by Eke |j] ■ 

The net effect is illustrated in figure 4, which shows 
the shift of the center of mass of the exciton transmission 
peak as a function of 7 for both theory and experiment. 
For wthz < ^12) we see a red shift which saturates and 
reverses with increasing intensity, showing a roughly lin- 
ear dependence on THz intensity at high fields (Fig. 4a). 
At low fields, A dominates, resulting in a net red shift. 
As A saturates with increasing field, Eke begins to dom- 
inate, eventually overwhelming the red shift entirely and 
resulting in a net blue shift. For wthz > ^12, we observe 
only an increasing blue-shift with increasing -Ethz (Fig. 
4b). Here, A and Eke cooperate to create a blue-shift. 

In conclusion, we have observed the experimental sig- 
nature of DFKE through its interplay with the AC Stark 
Effect. Our theory and observations are in remark- 
able agreement, and show the resulting frequency- and 
intensity-dependent shift of the excitonic resonance in a 
MQW under intense THz irradiation. 
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FIG. 1. A schematic map of electro-optic phase space as 
a function of amplitude and frequency of the perturbing field, 
showing DFKE regime "between" MPA and FKE. 




FIG. 2. (a) DFKE in an ideal 2D system. In a strong 
THz field (dashed curve, solid curve is zero-field), the "edge" 
shifts to higher frequency and sub-gap absorption increases, 
(b) The AC Stark Effect: a strong field applied at frequency 
uj ~ UJ12 causes 0^12 to shift. For u < (>)wi2, the transition 
shifts to w* 2 > (<Vi2- A = (u>\2 ~ uj^)/^ is the shift of 
|1) with respect to distant energy levels, (c) Schematic of 
experiment. 
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FIG. 3. Experimental transmission of MQW near elhl 
exciton with (a) Ticuthz = 2.5meV at 7thz = 0, 1, 2, 4, 12 (arb. 
units), (b) Tiwthz = 14meV at 7thz = 0, 1, 2, 4, 7 (arb. units). 
Arrows connect calculated centers of experimental peaks and 
point in direction of increasing 7thz- 
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FIG. 5. Calculated absorption of MQW with; (a) 
/iuthz = 2.5meV < hu)\2 from top down at 7 = 0, 0.03, 0.06, 
0.1, 0.5, 1.5. (b) ?iwthz = 14meV > huji2 from top down at 
7 = 0, 0.01, 0.02, 0.04, 0.08, 0.15. 
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FIG. 4. "Center of mass" of measured (squares) and cal- 
culated (triangles) exciton transmission peak vs. 7 for (a) 
Tllothz ~ 2.5meV < huii2 (b) and Swthz = 10.5meV > huji2. 
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